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ABSTRACT 

The clear characteristic timescale picked out by the low frequency quasi-periodic oscil- 
lations (QPOs) seen in many black hole and neutron star binaries has the potential to provide 
a very powerful diagnostic of the inner regions of the accretion flow. However, this potential 
cannot be realised without a quantitative model for the QPO. We have recently shown that 
the same truncated disc/hot inner flow geometry which is used to interpret the spectral transi- 
tions can also directly produce the QPO from Lense-Thirring (vertical) precession of the hot 
inner flow. This correctly predicts both the frequency and spectrum of the QPO, and the tight 
correlation of these properties with the total spectrum of the source via a changing truncation 
radius between the disc and hot flow. This model predicts a unique iron line signature as a 
vertically tilted flow illuminates different azimuths of the disc as it precesses. The iron line 
arising from this rotating illumination is blue shifted when the flow irradiates the approaching 
region of the spinning disc and red shifted when the flow irradiates the receding region of the 
disc. This gives rise to a characteristic rocking of the iron line on the QPO frequency which is 
a necessary (and probably sufficient) test of a Lense-Thirring origin. This is also an indepen- 
dent test of disc truncation models for the low/hard state, as vertical precession cannot occur 
if there is a disc in the midplane. 

We show that it may be possible to observe this effect using archival data from the Rossi 
X-ray timing explorer (RXTE) or XMM Newton. However, a clean test requires a combination 
of moderate resolution and good statistics, such as would be available from a long XMM- 
Newton observation or with data from the proposed ESA mission LOFT. 

Key words: X-rays: binaries - accretion, accretion discs 
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1 INTRODUCTION 

Low frequency quasi-periodic oscillations (QPOs) are commonly 
observed in the X-ray flux of both neutron star and black hole bi- 
naries (NSBs and BHBs respectively; collectively X-ray binaries 
XRBs). They are most clearly observed as strong, coherent features 
in the power spectral density (PSD) which are Lorentzian in shape 
and so can be described by amplitude (i.e. fractional rms variabil- 
ity), centroid frequency (/qpo) and width (A/). These properties 
are observed to be tightly correlated with the spectral properties of 
the source which vary dramatically as the source rises from quies- 
cence to outburst before falling once more into quiescence (see e.g. 
van der Klis 2006; Belloni 2010). 

The physical processes behind this spectral evolution are com- 
paratively well understood. The spectral energy distribution (SED) 
consists of three main components: a quasi-thermal disc, a power 
law (with high and low energy cut-offs) and a reflection spectrum. 
When the source flux is low (low/hard state), the power law is hard 
(photon index T ~ 1.7) and dominates the SED. As the source flux 
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increases, the power law softens and weakens while the disc and 
reflection spectra increase in luminosity (intermediate state). Even- 
tually at the peak of the outburst, the disc completely dominates 
(high/soft state), although sometimes there is also a strong high en- 
ergy tail (very high state). The disc spectrum is well explained by 
a standard thin disc (Shakura & Sunyaev 1973) and the power law 
can be reproduced by Compton up-scattering of disc photons by 
hot electrons in an optically thin corona. A fraction of the luminos- 
ity emitted from the corona will then reflect off the disc to give a 
reflection spectrum with the most obvious features being a strong 
iron K a line and a ~ 30 keV hump (see e.g. Fabian et al 2000). 

The truncated disc model, in which the thin disc only extends 
down to some radius r , can naturally explain the evolution of the 
SED (Esin et al 1997; Done, Gierlinski & Kubota 2007). Interior 
to r a is a large scale height, optically thin accretion flow (hereafter, 
the flow) which acts as the Comptonising corona. As the source 
flux increases, the truncation radius moves inwards, thus increas- 
ing the flux of disc photons incident on the flow and softening the 
power law emission while simultaneously decreasing all character- 
istic timescales associated with r . In this picture, r moves from 
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~ 60 — 6 (in units of R g = GM/c?) during the rise to outburst 
and back out again during the fall back to quiescence. 

QPOs are observed in the PSD during both the rise and the 
fall. During the rise, the PSD displays a QPO (with harmonics) 
superimposed on a broad band noise of variability. The broad band 
noise can be roughly characterised by two zero-centred Lorentzians 
with widths /(, and fh (e.g. Belloni, Psaltis & van der Klis 2002). 
The QPO frequency moves from ~ 0.1 — 10 Hz as the source 
flux increases and is correlated with rises in both ft, and fh (Wij- 
nands & van der Klis 1999; Psaltis, Belloni & van der Klis 1999). 
For BHBs, this is commonly classified as the type-C QPO. Even- 
tually the broad band noise disappears and the PSD is dominated 
by a type-B QPO which has fqpo ~ 6 Hz. Before the source 
completely transitions into the high/soft state, type-A QPOs are 
observed which are much broader and weaker features, typically 
centred at /qpo ~ 8 Hz. During the fall, the same is observed 
in reverse (see Casella, Belloni & Stella 2005 for more details of 
the A,B,C classification system and Belloni 2010 for a review of 
hysteresis behaviour). Because these three types of QPO are not 
observed simultaneously (even though the transition between type- 
A and type-B QPOs can be very rapid) and they occupy a similar 
frequency range, it is possible that they are driven by three different 
variants of the same underlying physical process. 

NSBs display a similar phenomenology of QPO types and 
spectral transitions (although the nomenclature is very different; 
see e.g. van der Klis 2005 for details). There is very strong evidence 
that the QPOs in both classes of object are produced by the same 
process. The correlations between fqpo, ft and fh hold, with the 
same gradient, for both NSBs and BHBs (Wijnands & van der Klis 
1999; Klein- Wolt & van der Klis 2008) with the only difference in 
frequency being entirely consistent with mass scaling (Ingram & 
Done 2011). 

The physical process responsible for driving the QPO remains 
poorly understood. There are many QPO mechanisms suggested 
in the literature (e.g. Fragile, Mathews & Wilson 2001; Wagoner 
et al 2001; Titarchuk & Oscherovich 1999; Cabanac et al 2010; 
Tagger & Pellat 1999; O'Neill et al 2011; Kato 2008; Wang et 
al 2012). However, such a rich phenomenology means that these 
models are rarely able to explain all of the known QPO proper- 
ties simultaneously. In Ingram, Done & Fragile (2009), we sug- 
gested perhaps the most promising QPO model to date. Based on 
the model of Stella & Vietri (1998), we considered the QPO to 
result from Lense-Thirring precession. This is a relativistic effect 
which occurs because a spinning compact object drags spacetime 
as it rotates. The orbit of a particle which is outside the plane of 
black hole spin will therefore undergo precession because the start- 
ing point of the orbit rotates around the compact object. Stella & 
Vietri (1998) and Stella, Vietri & Morsink (1999) showed that the 
Lense-Thirring precession frequency of a test mass at the truncation 
radius is broadly consistent with the QPO frequency. Schnittman 
(2005) and Schnittman, Homan & Miller (2006) developed this 
into a fully relativistic description of a misaligned ring, showing 
that its direct emission and iron line signature should be modu- 
lated on the precession frequency, which could be somewhat higher 
than observed. However, the real problem with these models is that 
the energy spectrum of the QPO is dominated by the Comptonised 
emission (Sobolewska & Zycki 2006; Rodriguez et al 2004), re- 
quiring that the QPO mechanism predominantly modulates the hot 
flow rather than the disc (although the variability could be produced 
elsewhere before propagating into the flow; Wilkinson 2011). We 
consider instead a global precession of the entire hot flow, which 
naturally explains the QPO spectrum. Such global precession has 



been seen in recent numerical simulations (Fragile et al 2007, Frag- 
ile 2009). We show in Ingram, Done & Fragile (2009) that the pre- 
dicted frequency range is completely consistent with the type-C 
QPO in BHBs and also in NSBs (Ingram & Done 2010). 

There are other more subtle properties that are naturally pre- 
dicted by the precessing flow model. Heil, Vaughan & Uttley 
(201 1) show that the QPO frequency is linearly related to the source 
flux on short timescales (~ 3s). We show in Ingram & Done (201 1) 
that propagating fluctuations in mass accretion rate which give rise 
to the broad band noise (e.g. Lybarskii 1997; Arevalo & Uttley 
2006) will affect the moment of inertia of the flow leading the pre- 
cession frequency to fluctuate. The linear relation with flux then 
occurs because both the flux and the precession frequency depend 
on mass accretion rate. Although it is very encouraging that this 
property is predicted by the model, we still do not have unambigu- 
ous proof that the flow precesses - a QPO produced from any mode 
of the hot flow will also couple to fluctuations propagating through 
the hot flow, and should give an /qpo -flux relation. 

The interpretation of the QPO as vertical precession requires 
a truncated disc as otherwise the flow could not cross the equatorial 
plane. The issue of whether or not the disc truncates is still some- 
what controversial. The line clearly depends on the spectral state, 
with a very small narrow line seen in the dimmest low/hard states 
(e.g. Tomsick et al 2009), and a very broad line when the source 
is very close to the transition to the soft state (e.g. Hiemstra et al 
2011; a hard intermediate state just after the transition from the 
soft state). However, in the brighter low/hard states, Nowak et al 
(2011) show that the broad iron line in Cyg X-l can be variously 
interpreted as implying a disc anywhere from 6 — 32R g (for their 
Obs 4) depending whether the continuum is thermal Comptonisa- 
tion, non-thermal Comptonisation, multiple Compton components 
or includes a jet contribution. Fabian et al (2012) show another de- 
convolution of a similarly shaped spectrum from Cyg X-l, where 
the spectrum below 10 keV is dominated by highly ionised, highly 
smeared reflection, with a very small inner radius of ~ 1.3i? 9 and 
a very steep emissivity profile (a.k.a. the lightbending model). We 
note that this lightbending geometry is inconsistent with the in- 
dependent requirement on the un-truncated disc geometry that the 
source is beamed away from the disc in order to produce an intrin- 
sically hard spectrum (Malzac, Beloborodov & Poutanen 2001). 

The issue is clearly still very controversial, though we note 
that the rapid spectral variability can only currently be explained 
with an inhomogeneous Comptonisation continuum model (Kotov 
et al 2001; Arevalo & Uttley 2006), where the line profile is consis- 
tent with a truncated disc (Makishima et al 2008). Here we simply 
assume the truncated disc geometry, and use this to propose a dis- 
tinctive test of a vertical precession origin of the QPO. As a tilted 
flow precesses, the illumination pattern on the disc rotates. The re- 
sulting iron line is boosted and blue shifted at a time when the flow 
illuminates the approaching side of the disc, and red shifted when 
the flow illuminates the receding side of the disc. Since this peri- 
odic rocking of the iron line is a requirement of the Lense-Thirring 
QPO model, this also offers a potentially unambiguous test of disc 
truncation. Our geometry differs from the Schnittman et al (2006) 
model, where a precessing inner disc ring producing the iron line 
and continuum. Instead, we have a hot inner flow replacing the in- 
ner disc to produce the continuum, and precession of the entire hot 
flow produces a rotating illumination pattern which excites the iron 
line from the outer thin disc. 

The paper is ordered as follows. In section [2] we define the 
accretion geometry assumed for the model. In section [3] we will 
calculate the implications of our assumed geometry on a very sim- 
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pie toy spectral model. In section [4] we will take this further by 
introducing a reasonable spectral model before analysing the like- 
lihood of observing this effect in section [5] 



2 MODEL GEOMETRY 

In this section, we outline the geometry used for our QPO model. 
We assume that the spin axis of the compact object is misaligned 
with that of the binary system as may be expected from supernova 
kicks (Fragos et al 2010). Due to frame dragging, the orbit of an 
accreting particle from the binary partner will precess around the 
spin axis of the compact object. The effect of frame dragging on an 
entire accretion flow depends on the dynamics of the flow. A thin 
accretion disc being fed by a binary partner out of the spin plane 
of the compact object will form a Bardeen Petterson configuration 
(Bardeen & Petterson 1975) where the outer regions align with the 
binary partner and the inner regions align with the spin of the com- 
pact object, with a transition between the two regimes at rpp- The 
value of tbp is not well known, with analytical estimates ranging 
from ~ 10 - 400 R 9 (see e.g. Bardeen & Petterson 1975; Pa- 
paloizou & Pringle 1983; Fragile, Mathews & Wilson 2001). In the 
thin disc regime, warps caused by the misaligned black hole prop- 
agate in a viscous manner. This means that the time scale on which 
a warp is communicated is much longer than the precession period 
and therefore a steady configuration forms. In contrast, warps in 
a large scale height accretion flow are communicated by bending 
waves (see e.g. Lubow, Ogilvie & Pringle 2002; Fragile et al 2007) 
which propagate on approximately the sound crossing timescale 
which is shorter than the precession period. For this reason, the hot 
flow can precess as a solid body with the precession period given 
by a surface density weighted average of the point particle preces- 
sion period at each radius (Liu & Melia 2002), while a cool disc 
forms a stable warped configuration. This solid body precession of 
a hot flow has been seen explicitally in recent numerical simula- 
tions (Fragile et al 2007) for the special case of a large scale height 
flow which we consider here. 

The key aspect is that the flow angular momentum has to be 
misaligned with the black hole spin. Yet the outer thin disc will 
warp into alignment with the black hole at rpp- Since this radius 
is poorly known, there are two possible scenarios. Firstly rpp may 
be small enough for the outer thin disc to still be aligned with the 
binary partner at the truncation radius. In this case, the hot flow is 
misaligned with the black hole spin by the intrinsic misalignment 
of the binary system which will naturally lead to solid body preces- 
sion of the entire flow. Secondly, if rpp is large, the disc and hence 
hot flow are intrinsically aligned with the black hole spin. However, 
precession may be possible. The flow has a large scale height, so is 
sub-Keplarian. At the truncation radius it overlaps with the Keplar- 
ian disc, so this overlap layer is probably Kelvin-Helmholtz unsta- 
ble, producing turbulence. Clumps forming from random density 
fluctuations in regions high above the midplane could temporarily 
mis-align the flow leading to intermittent precession. This predicted 
intermitency has the advantage of naturally explaining the observed 
random jumps in QPO phase (Miller & Homan 2005; Lachowicz 
& Done 2010). 

Here we assume the first geometry i.e. assume that rpp is very 
small. However, the effect of rotating illumination on the iron line 
is qualitatively the same in the second geometry, differing only in 
the details. In the next section, we outline the geometry used. We 
work under the assumption that the central object is a black hole, 
but the geometry is valid for neutron stars also. 



2.1 Disc 

The geometry we consider for the two component accretion flow 
is illustrated in Figure Q] We assume that the disc has angular mo- 
mentum vector set by the binary system, J_ BS , and that this is mis- 
aligned with the spin axis of the black hole (the z-axis) by an angle 
P. The flow angular momentum vector, J_f iom , precesses around 
the z-axis with phase given by the precession angle, 7. The plane 
of the disc is the plane orthogonal to j_ BS , while the plane of the 
flow is orthogonal to j_n w ^ n m i s coordinate system, the binary 
partner will orbit in the 'disc' plane. The observer's position is de- 
scribed by an inclination angle, 6i, and a viewer azimuth, which 
can take the range of values ^ 8i ^ n/2 and ^ <j>i ^ 27r. Here, 
6i is defined with respect to the binary (i.e. the disc) angular mo- 
mentum vector and (f>i is defined with respect to the x-axis. 

The flow then precesses around a circle centred on the black 
hole spin axis, from being aligned with the disc when 7 = 0, to 
being misaligned by angle 2/3 with respect to the disc when 7 = tt. 
We can define a vector r d which points from the black hole to any 
point on the disc plane. If the top of the flow is its brightest part, 
the region of the disc most strongly illuminated by the flow for a 
given 7 is where the angle between r d and J_f low is smallest. The 
smallest this angle can ever be is for r d = e when 7 = tt; i.e. 
this is the most that the flow angular momentum vector ever aligns 
with any azimuth of the disc plane, e therefore defines the azimuth 
of the disc which sees the maximum illumination from the flow. 
Material in the disc is spinning rapidly and, because precession is 
prograde, this orbital motion is anti-clockwise for our geometry. 
The viewer azimuth <j>i therefore specifies the direction with respect 
to the viewer in which disc material in the maximally illuminated 
region (i.e. on the e axis) is moving. For <f>i = 0, the receding part of 
the disc is most strongly illuminated as the flow precesses. Instead, 
for cj>i = 7r/2 the maximum illumination is on the patch directly 
in front of the black hole. For (f>t — tt the maximum illumination 
is on the approaching side of the disc, while for 3n/2 it it for the 
patch directly behind the black hole. We assume that the disc is 
razor thin and flat (i.e. no flaring). The mathematical definitions for 
the geometry we use are outlined in Appendix [A] 



2.2 Flow 

Unlike the disc, the flow has volume and scale height, so we must 
make some assumptions about its shape. We assume that it is a 
crushed sphere; i.e. viewed from above it is circular but it has 
an elliptical cross section as illustrated in Figure [2] The semi- 
major axis of the ellipse is r and the semi-minor axis is h . We 
choose to parametrise this by defining a scaleheight, h/r, such that 
h = (h/r)r a . Figure [2] also shows that we set an inner radius, n, 
such that the core of the quasi-spherical flow is missing. This is to 
incorporate a flavour of the numerical simulations which show that 
shocks (at the bending wave radius) can truncate the inner region 
of the hot flow (Fragile et al 2007). Any point on the flow surface 
is then a distance r away from the black hole, where r is a function 
of the angle 9f. We assume that each radius of the surface radiates 
the gravitational potential energy released at that radius (i.e. we 
use a surface rather than a volume emissivity). This gives a simple 
analytic model where the central parts of the flow (outside of rj) 
are brighter than the outer parts, but that these bright regions are 
near the poles which gives a reasonable reflection fraction, (Q/2tt) 
while also giving a reasonable precession frequency (set by r;, r D , 
Al, the surface density profile which we assume to be constant, and 
a», where a* is the dimensionless spin parameter: equation 1 in In- 
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Figure 1. Schematic diagram illustrating the coordinate system we are con- 
sidering. The black hole is at the origin and the black hole angular momen- 
tum vector is aligned with the z-axis. J s g is the (unit) angular momentum 
vector of the binary system (and the disc), misaligned with the z-axis by 
an angle ft. £ then completes a right handed Cartesian coordinate system 
{x, e, J_ B g} such that the disc plane is described by the plane J_ BS = 0. 
J_fl ow is the angular momentum vector of the flow and we see its orienta- 
tion precesses around the blue dotted ring, its phase described by the pre- 
cession angle 7. The flow, shown in (translucent) blue, is then described by 
the plane orthogonal to J_fi ow - The observer's position is described by 6>; 
and (pi. 




Figure 2. Schematic diagram illustrating the cross section of the flow. See 
text for details. 

gram, Done & Fragile 2009). Note that, eventhough this is a simpli- 
fied prescription, the most influential aspect of the flow geometry is 
where the brightest region lies. In nearly all imaginable geometries, 
this point lies at the pole of the flow (as it does for our geometry). 
Thus our mathematically convenient assumptions for flow geome- 
try should provide us with results not materially different from a 
far more difficult calculation assuming a geometry identical to the 
Fragile et al (2007) simulation. More details of the flow geometry 
are presented in Appendix lAl 

Fundamentally, the precession frequency modulates the con- 
tinuum as the pole moves in and out of sight. The QPO maximum 
occurs when the pole faces the observer and the minimum when it 
faces away. Thus the region of the disc preferentially illuminated is 
in front of the black hole (from the point of view of the observer) 
at the QPO maximum and behind for a QPO minimum. Because 
precession is prograde, this means that the flow illuminates the ap- 
proaching disc material during the rise to a QPO maximum (be- 
cause the pole has to first move towards us in order to face us) and 
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Figure 3. Disc irradiation by the flow as seen by a viewer with 0i = 60° 
and (pi = 0° (left) or (pi = 90° (right). The flow is shown in grey with 
black gridlines for clarity. The truncation radius is r = 60. The luminosity 
incident on the disc is grouped into 8 bins with black, red, green, blue, cyan, 
magenta, yellow and orange representing the dimmest to brightest patches 
on the disc. The solid black line in the top picture of each plot indicates the 
black hole spin axis. Flow precession causes the characteristic illumination 
pattern to rotate around the disc. 

the receding material on the fall to a QPO minimum. Below we 
calculate the self-consistent illumination pattern for the disc as a 
function of QPO phase for our assumed geometry. 



3 IMPLICATIONS OF A PRECESSING FLOW 
3.1 Disc irradiation 

Each flow surface element will radiate a luminosity dL over a semi- 
sphere (because the element radiates away from the black hole). A 
disc surface element with area dAd will intercept some fraction 
of this luminosity. This fraction can be calculated self-consistently 
from the projected area of the disc element. The disc element will 
not intercept any of the luminosity from the flow element if it makes 
an angle greater than 7r/2 with a vector which is orthogonal to the 
flow element and points away from the black hole (i.e. if it is not in 
the unit semi-sphere of the flow element). Also, as observers with 
8i ^ 90°, we only observe reflected photons which are intercepted 
by the top of the disc. 

The total incident luminosity on the disc surface element is 
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Figure 4. The iron line profile as seen by a viewer with 8i = 60° and 
4>i = 0° (left) or <j>i = 90° (right). The rest frame iron line profile is 
assumed to be a S— function at 6.4 keV and the truncation radius is r a = 60 
as in Figure [3] Different colours represent different snapshots in time with 
black, red, green, blue and cyan representing the top to bottom snapshots 
pictured in Figure[5] The rotation of the illumination pattern causes the iron 
line profile to rock from red to blue shift. 




5 5.5 6 6.5 7 5 5.5 6 6.5 7 

Energy (keV) Energy (keV) 

Figure 5. The iron line profile as seen by a viewer with 9i = 60° and 
4>i = 0° (left) or 4>i = 90° (right). The rest frame iron line profile is still 
assumed to be a S— function at 6.4 keV but the truncation radius is now 
r = 10. The different colours represent the same snapshots in time as in 
Figure \4\ We see the motion of the iron line is different here compared to 
Figure|4] Due to stronger Doppler and relativistic boosting in the inner disc, 
the red wing never dominates in the E < 6.4 keV region. 



calculated by integrating over the surface of the entire flow. We do 
this calculation for every disc surface element over a full range of 
precession angles (0 ^ 7 < 2n) in order to build a picture of disc 
irradiation as a function of precession angle (and therefore time). 
The details of this calculation are presented in Appendix [B] For 
simplicity, we use a Euclidean metric i.e. assume that light travels 
in straight lines. This should be a fairly reasonable approximation 
because we assume a fairly large value of Ti throughout the paper 
(following Dexter & Fragile 2011; Ingram & Done 2012; Fragile 
2009) and so lightbending is not very significant (e.g. Fabian et al 
1989). 

Throughout the paper, we will use the values r , = 7, /3 = 15° 
and h/r = 0.9 (we discuss our reasoning for these fiducial values 
in section [3~3b . Figure [3] shows the resulting illumination pattern 
with r a = 60, with snapshots taken at five different values of pre- 
cession angle 7 for an inclination angle of 6i — 60°. The left hand 
plot shows the pattern as seen by an observer at cj>i — 0°, whereas 
the right hand plot shows this for 4>i = 90°. The luminosity is 
grouped into bins of equal logarithmic size with black, red, green, 
blue, cyan, magenta, yellow and orange representing the dimmest 
to brightest bins respectively. The flow is shown in grey with black 
gridlines included for clarity. In the top picture of each plot, we also 
include a straight black line to illustrate the orientation of the black 
hole spin axis. This is misaligned with j_ BS by /3 = 15° but, as 
Figure [3] demonstrates, the apparent misalignment between these 
two vectors depends on the viewing position. We clearly see the 
flow precess, with the pole of the flow moving in a circle around 
the black hole spin axis. As it does, the brightest part of the disc 
is always the region closest to the pole of the flow meaning that 
it rotates around the disc. Because of our asymmetric geometry, 
the flow starts off aligned with the disc, is misaligned by 2/3 when 
7 = 7r before aligning again for 7 = 2n. For cj>i = 0, the maximum 
misalignment (giving the maximum illumination of the disc) is on 
the right hand (receding) side of the disc, while for 0, = 90 it is 
directly in front of the black hole, but in both cases the illumination 
pattern rotates. In the next section, we will discuss how this will 
affect the observed iron Kq, line. 

3.2 Effect on the iron K Q line profile 

When the flow emission irradiates the disc, bound atoms in the disc 
will fluoresce to produce emission lines, the most prominent being 
the iron K a line at ~ 6.4 keV (George & Fabian 1991; Matt, Per- 
ola & Piro 1991). However, this line is in the rest frame of the disc 
which is rotating rapidly meaning that a non face-on observer will 
see some regions of the disc moving towards them and others reced- 
ing. Doppler shifts mean that emission from the approaching side 
is blue shifted while that on the receding side is red shifted. Also, 
length contraction along the line of motion beams the emission in 
that direction. Thus the blue shifted emission from the approach- 
ing side is also boosted in comparison to the red shifted emission, 
leading to a broadened and skewed iron line. An additional energy 
shift is provided by time dilation and also gravitational redshift 
which combine to broaden the line even further (Fabian et al 1989; 
2000). Figure [3] clearly shows that, according to this model, the 
disc irradiation pattern rotates around the disc meaning that some- 
times the brightest region of the disc is receding (e.g. the </>i = 0°, 
7 = 47r/5 scenario in Figure[3), and sometimes the brightest region 
is approaching (e.g. the (j>i — 90°, 7 = 27r/5 scenario in Figure 
[3J. Therefore, as the flow precesses, the iron line will periodically 
rock between red and blue shift. In this example, the material in 
the disc and the irradiation pattern are both rotating anti-clockwise. 
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In general, they they could both be moving clockwise but the re- 
sulting pattern is the same (maximum blueshift, QPO maximum, 
maximum redshift, QPO minimum). Lense-Thirring precession is 
prograde, so the disc and flow will never be rotating in opposite 
directions, making this periodic shifting of the iron line profile a 
unique prediction of the model. 

We use the illumination pattern on each surface element of the 
disc to set the amount of intrinsic iron line emission. We assume 
that this is a 8— function at Eq = 6.4keV and then use the radius 
and azimuth of the surface element of the disc and the inclination of 
the observer to calculate the shifted line emission (see Appendix). 

Figure [4] shows the iron line profile at five snapshots of time 
with black, red, green, blue and cyan lines corresponding to 7 = 0, 
2n/5, 4/7r/5, 6n/5 and 8n/5 respectively. We use the same param- 
eters as for Figure [3] The details of this calculation are presented 
in Appendix [D] For simplicity, we do not include light bending but 
this should not be a large effect for the comparatively large radii 
we consider. The left hand plot is for <j>i — 0°, r — 60 (i.e. cor- 
responding to the left plot of Figure O and we see that the iron 
line does indeed rock between red and blue shift as the illumina- 
tion pattern rotates. Note that, for these parameters, the 2nd and 
5th snapshots have an identical iron line profile, as do the 3rd and 
4th snapshots. The right hand plot is for <j>i — 90°, r a = 60 (i.e. 
corresponding to the right hand plot in Figure [3}. We see that the 
periodic rocking has a different phase and the peak flux of the blue 
wing is much larger. This is because, for the cj>i = 0° case, the ap- 
proaching side of the disc is never the brightest part, whereas this 
does happen for the <j>i = 90° case. This movement of the iron 
line is obviously a very distinctive model prediction and so could 
provide a detectable, unambiguous signature of a vertically tilted, 
prograde precessing flow i.e. a clean test of a Lense-Thirring origin 
of the QPO. 

Figure [5] shows the same thing but now r = 10. We see that 
Doppler (and relativistic) boosting of the blue wing is now such a 
large effect that the red wing never dominates even when the flow is 
preferentially illuminating the receding material. As such, the mo- 
tion of the iron line is different. Crucially, although the exact shape 
of the iron line depends on the illumination pattern and thus the 
details of the assumed flow geometry, this dependence on trunca- 
tion radius is really quite robust to changes in flow geometry. The 
differences between Figures[4]and[5]are driven primarily by the dif- 
ference in disc angular velocity (i.e. the position of the truncation 
radius) and not the details of the modelling. Thus this effect could 
provide a robust diagnostic for the accretion flow geometry. 



3.3 Modulation of the continuum 

As the flow precesses, the luminosity seen by the observer will 
change periodically giving rise to a strong QPO (with the quasi- 
periodicity provided by frequency jitter among other processes; 
Ingram & Done 2012; Heil, Vaughan & Uttley 2011; Lachowicz 
& Done 2010). This is because the total surface area of the flow 
viewed by the observer changes and, also, some regions of the flow 
are brighter than others meaning that a trough in the light curve 
would typically occur when the brightest regions of the flow (i.e. 
the poles) are hidden. The calculation for this process is similar to 
that performed in section [3~T1 Each flow surface element emits a 
luminosity dL. The observer at 8i, (pi will see no luminosity from 
this surface element if they are not within the unit semi-sphere of 
the element, and we also remove luminosity from lines of sight 
which are obstructed by the disc. We can then integrate over every 



flow element to calculate the observed luminosity as a function of 
precession angle. 

The blue lines in Figure [6] show the observed luminosity ex- 
pressed as a fraction of the total luminosity, L to t, plotted against 
precession angle. We use the fiducial parameters n = 7, ft = 15° 
and h/r — 0.9 and consider the r = 60 example. The solid line 
is for (j>i — 0° and the dashed line represents (pi — 90°. As ex- 
pected, the observed luminosity varies with precession angle and 
the phase depends on <j>i. The fractional rms is 8.4% and 4.2% for 
(pi — 0° and (pi = 90° respectively. These values are lower than 
the observed QPO rms values of ~ 10 — 15%. However, the pre- 
dicted values would be higher if we were to consider that the flow 
is fed by disc photons, the flux of which incident on the flow will 
change periodically as the flow precesses. We ignore this process 
here because it will affect the direct and reflected emission equally 
and so will not contribute to the rocking iron line effect. 

For the green line, we plot the total luminosity incident on the 
disc (which determines the iron line / reflected flux) as a function of 
precession angle. Because the disc is flat, this does not depend on 
(pi. This effectively tracks the misalignment between flow and disc 
with the minimum reflection occurring when the flow is aligned 
(7 = 0) and the maximum when the flow is misaligned by 2/3 
(7 = 7r). Hence the direct and intercepted emission are generally 
out of phase. The black lines show the reflection fraction (inter- 
cepted/direct) with the solid and dashed lines representing (pi = 
and 90° respectively. This corresponds to the solid angle of the 
disc, and the time averaged ratio for (j>i — 0° is Q,/2n = 0.263, 
and with Q/2n = 0.238 for the (j>i = 90° case. These values are 
fairly representative of those observed for the low/hard state (e.g. 
Gierlinski et al 1999; Zycki, Done & Smith 1998; Gilfanov 2010). 

Note that large value of h/r gives a reasonable reflection frac- 
tion but under predicts the QPO rms. If we had considered, for ex- 
ample, an overlap region between disc and flow, disc flares or a 
small disc scale height, we could have achieved a reasonable re- 
flection fraction and the correct QPO rms (for this we would also 
need to consider the variation in disc seed photons) for a far lower 
value of h/r. However, these effects are all very difficult to model 
and our assumed geometry should not significantly affect the final 
results. Thus we choose the fiducial parameter values to give rea- 
sonable results for a simplified geometry. 



4 SPECTRAL MODELLING 

We now use a full reflected spectrum rather than just a line, and 
recalculate the effect of the rotating disc illumination pattern and 
varying effective area of the flow for this more realistic scenario. 
We consider the same two values of truncation radius as those con- 
sidered previously, r = 60 and r — 10. These values corre- 
spond to precession frequencies of fQpo(r = 60) = 0.145 Hz 
and fQpo(f — 10) = 5.36 Hz for the fiducial parameters, a 
spin of a, = 0.5 and a mass of M = IOMq (i.e. 2.9 and 107.1 
a* (Mq /M) Hz). They also correspond to different spectral states, 
with r = 60 giving rise to a low/hard state (LHS) spectrum and 
r — 10 leading to a soft intermediate state (SIMS) spectrum. The 
QPO in the LHS spectrum will be of type-C whereas it will be of 
type-B for the SIMS spectrum. 

4.1 Method 

For both the LHS and SIMS spectra, we include quasi-thermal disc 
emission, Comptonised flow emission and a reflection spectrum. 
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Figure 6. Emission as seen by a viewer at Oi = 60° and <f>i = 0° (solid 
lines) or </>; = 90° (dashed lines). The blue line represents emission di- 
rectly observed from the flow. We see that precession of the flow introduces 
a strong periodicity. The green line represents the total luminosity inter- 
cepted by the disc. This also has a periodicity because the misalignment 
between disc and flow changes as the precession angle, 7, evolves. It does 
not, however, depend on the position of the observer. The black line is the 
ratio between direct and reflected (intercepted) light, (Q/2tt). 



We use XSPECvl2 (Arnaud 1996) throughout. We describe the disc 
with DISKBB (Mitsuda et al 1984), and for simplicity we assume 
that this spectrum is constant. This is not strictly true. Figure [3] 
shows that the inner disc is periodically obstructed by the flow, giv- 
ing a small periodicity in the hottest part of the disc emission. Also, 
the non-reflected photons which illuminate the disc will thermalise 
and add to the intrinsic disc emission, and this additional thermal 
emission will vary in intensity, being stronger when the flow is at 
its maximum misalignment angle to the disc, and weakest when 
the flow is aligned with the disc. This additional thermal emission 
is also periodically redshifted/blueshifted in the same way as the 
line. However, these effects should be small as they are diluted by 
the much larger constant flux from the disc. We will investigate 
this in a future paper, as evidence for this may have been observed 
(Wilkinson 201 1). However, here we are interested in the iron line 
region and so ignore this potential contribution to the QPO in the 
disc spectrum. 

For the flow we assume that every element emits the same 
spectrum, meaning that the periodicity is in the normalisation of the 
flow spectrum. We describe the spectrum by the Comptonisation 
model NTHCOMP (Zdziarski, Johnson & Magdziarz 1996; Zycki, 
Done & Smith 1999) which produces a power law spectrum with 
high and low energy cut-offs governed by the electron temperature 
and disc photon temperature (kTbb tied to the disc temperature) 
respectively. We fix the normalisation of this by the angle averaged 
flux from the flow (Ltot), to set the flux from each surface element 
of the flow. We then use the method described in section [3~3l to 
determine the modulation of the observed continuum, to calculate 
the factor by which to multiply the normalisation of NTHCOMP as 
a function of phase angle. 

We use the method described in section [3~2l to calculate the 
illuminating flux from the flow at each surface element in the disc, 
and use this to set the normalisation of the illuminating NTHCOMP 
model. We describe the shape of the resulting reflected emission 
by RFXCONV (Ross & Fabian 2005; Done & Gierlinski 2006; 
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Table 1. Summary of the parameters used for both the LHS and SIMS spec- 
tral models. 

Magdziarz & Zdziarski 1995; Kolehmainen, Done & Diaz Trigo 
2011). This is similar in form to the IREFLECT model in XSPEC 
but replaces the very approximate ionisation balance incorporated 
in this model with the much better Ross & Fabian (2005) calcula- 
tions. This outputs a partially ionised (parametrised by log 10 £ ) re- 
flection spectrum, including the self consistent emission lines, for 
a general illuminating spectrum. We fix the inclination angle of the 
reflector at Qi and abundances at solar. We calculate the reflected 
emission from this illuminating flux assuming £l/2n — 1. This is 
an underestimate as RFXCONV assumes that the disc is illuminated 
isotropically, whereas in our geometry the illumination is prefer- 
entially at grazing incidence. However, the amount of reflection is 
also set by the unknown details of the shape of the flow, so this ap- 
proximation is good enough to demonstrate the general behaviour 
of the model. 

The reflected emission from each surface element is shifted 
in energy depending on the radius and azimuth (see Appendix fDt. 
We sum the reflected emission from all the disc elements to derive 
the total reflected emission for each phase. This gives the correct 
relative normalisation of the continuum and reflected flux, and how 
this changes as a function of precession phase angle 7 for a given 
set of model (r , r^, /?, h/r, 9i, (f>i) and spectral (kT bb , T, log£, 
kT e ) parameters. 

4.2 Phase resolved spectra 

The parameters used for each state are shown in Table Q] We as- 
sume that kT bb , r and disc ionisation increase as the rise to out- 
burst continues whereas kT e decreases, as is commonly observed. 
The resultant time averaged LHS (r = 60) spectrum has a 2 — 10 
keV flux of ~ 0.3 Crab and Q/2n — 0.24 (or iron line equivalent 
with of 150 eV when fit by a diskline profile rather than a full re- 
flected spectrum). For the SIMS (r = 10) spectrum, the flux is 
~ 0.66 Crab and the reflection has Q/2n — 0.42 (iron line equiv- 
alent width of ~ 240 eV) with a much steeper continuum. These 
values are typical of those observed in the relevant states for fairly 
bright BHBs (e.g. GRS 1915+104 in its QPO state: Ueda et al 2010, 
and the intermediate state of GX339-4; Tamura et al 2012), justify- 
ing our choice of parameters. 

Figure [JJ shows the LHS spectrum as viewed from a position 
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Figure 7. LHS spectrum for five snapshots in time calculated using the 
model described in the text, using the parameters listed in Table[T] We use 
the same convention as for Figures[4]and[5]wifh black, red, green, blue and 
cyan representing the first to last snapshots. The top plot is a broadband 
spectrum with all of the components. The disc and Comptonisation compo- 
nents are both represented by dotted lines and the total spectrum as well as 
the reflection component are represented by solid lines. The bottom right 
plot zooms in on the intrinsic iron line and the bottom left plot zooms in on 
the iron line region of the total spectrum. We see that the motion of the iron 
line is still present but dilution from the continuum makes the effect much 
more subtle in the total spectrum. 



with <j>i = 90° and 9i = 60° at five different snapshots in time. 
We use the same convention as for Figures[4]and[5]with black, red, 
green, blue and cyan representing 7 = 0, 2n/5, Atv/5, 6ir/5 and 
8n/5 respectively. The top plot shows the total spectrum (upper 
solid lines) and its components, the constant disc (black dotted line 
just seen in the lower left hand corner of the plot), variable flow 
(dotted continuum lines just underneath the total spectra - the sym- 
metry means that the red dotted line is the same as the cyan, while 
the green is the same as the blue) and reflected spectra (lower solid 
lines). We clearly see the flow continuum oscillate while the reflec- 
tion spectrum rocks between red and blue shift, as well as changing 
in normalisation. The reflection spectrum is in phase with the con- 
tinuum in this example because <f>i — 90° (see Figure [6j but, in 
general, there is a phase difference between the two components. 
The lower left plot zooms in on the iron line region in the total spec- 
trum, while the lower right plot shows the changes in the reflected 
emission. We see that the reflected spectrum displays similar be- 
haviour to the corresponding 8— function (right hand plot of Figure 
[4^. The rocking movement in the underlying reflection spectrum is 
still visible in the total spectrum, though somewhat diluted by the 
changing continuum level. 

Figure [8] shows the same thing but for the SIMS. As for the 
S— function iron line profile in section fX2l we see that the major 
effect is now the strength and position of the blue wing rather than 
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Figure 8. SIMS spectrum for five snapshots in time calculated using the 
model described in the text, using the parameters listed in Table[T] We use 
the same conventions as for Figure [7] We see that, as for the 8— function 
calculation, the movement of the iron line is characteristically different for 
the SIMS compared with the LHS. 

a rocking motion from blue to red due to the much stronger Doppler 
and relativistic boosting in the inner disc. Non theless, there is still 
a clear periodic shift in the line shape with QPO phase, although 
the pronounced rocking of the iron line peak energy predicted for 
the LHS provides more of a 'smoking gun' for the Lense-Thirring 
model. 



5 OBSERVATIONAL PREDICTIONS 

In this section, we consider how this effect may be best observed. 
One potential method is to look at phase lags between differ- 
ent energy bands. We could define a red wing energy band (say 
5.4 — 6.4keV) and a blue wing energy band (say 6.4 — 7.4keV) and 
look for a phase lag between the two. However, Figures [7] and [8] 
show that, due to dilution from the periodically varying continuum, 
the energy shifting of the iron line is very subtle in the total spec- 
trum. This means that the phase lag between red and blue wings 
is very small (2 — 6 x 10 _2 7r) for our model and, consequently 
may be difficult to observe. Instead, we consider phase resolved 
spectroscopy. 

5.1 Phase binning 

The random phase jumps and varying period characteristic of QPO 
light curves make phase resolved spectroscopy difficult. Naively 
folding the light curve on the QPO period is not appropriate. It 
is, however, possible to isolate the maximum and minimum phase 
bins of the QPO by averaging over the brightest and faintest points 
in the light curve. Miller & Homan (2005) did this for two GRS 
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Figure 9. The 2-20 keV integrated flux of the LHS model with </>; = 90° 
and 9i = 60° plotted against precession angle. The dashed lines are flux 
thresholds. Intervals of the light curve above the top dashed line are con- 
sidered to be the QPO peak, intervals below the bottom dashed line are 
considered to be the trough. The rising section which will always follow 
a trough will have the bluest iron line profile. The falling section which 
always follows the peak will have the reddest iron line profile. 

1915+105 light curves, both containing a strong type-C QPO. This 
allowed them to compare the spectra corresponding to the QPO 
peak and trough. This analysis can be taken a step further because a 
rise will always follow a trough and a fall will always follow a peak. 
This simple phase binning can therefore provide four phase bins as 
opposed to two. Crucially, our model predicts that the maximum 
red shift always follows the QPO peak and the maximum blue shift 
always follows the QPO trough. This is because the pole of the flow 
(which is the brightest region) faces us, then is moving away from 
us, then faces away from us, then is moving towards us (before 
facing us again). Therefore the flow illuminates the observer, then 
the receding (red shift) part of the disc, then the region hidden to 
the observer then the approaching (blue shift) part of the disc. 

Figure |9] shows the 2-20keV light curve of our LHS model 
with 4>i = 90° and 9i — 60°. We define a peak as the bright- 
est 10% of the light curve and a trough as the faintest 10%. These 
thresholds are shown as dashed lines. We can therefore isolate the 
trough, the blue rise, the peak and the red fall. This flux selection 
means that the majority of the counts lie in the more interesting rise 
and fall sections as opposed to the peak and trough (unlike the flux 
selection of Miller & Homan who were interested in the peak and 
trough spectra). Figure[l0](top) shows the result of averaging spec- 
tra belonging to each of these four phase bins. The green line is the 
trough spectrum, the blue line is the rise spectrum, the black line is 
the peak spectrum and the red line is the fall spectrum. All are plot- 
ted as a ratio to a power law with photon index r = 1.6. We use 
this photon index rather than F — 1.7 because the reflection hump 
makes the total spectrum harder than the underlying Comptonisa- 
tion. As expected, the rise spectrum contains the most heavily blue 
shifted iron line and the fall spectrum contains the most heavily red 
shifted iron line. Because we tie the normalisation of the power law 
across the four spectra, we can see that the peak spectrum has the 
highest flux, the trough spectrum has the lowest and the rise and 
fall have comparable flux. 

In the bottom plot of Figure [Tol we plot the red fall spectrum 
subtracted from the blue rise spectrum. We use the absolute spec- 
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Figure 10. Top: Phase binned spectra calculated assuming r = 60, 4>i = 
90° and 0; = 60° plotted as a ratio to a power law with photon index T = 
1.6. These four phase bins are for the QPO minimum (green), rise (blue), 
maximum (black) and fall (red). As expected, the rise has the most heavily 
blue shifted iron line and the fall has the most heavily red shifted iron line. 
Bottom: The red fall spectrum subtracted from the blue rise spectrum. The 
solid line is for the r = 60 example shown in the top plot and the dashed 
line is for r D = 10. The shape of this difference spectrum is different for the 
two truncation radii. There is no negative section in the dashed line because 
strong Doppler and relativistic boosting in the inner disc prevents the red 
wing from dominating. 

trum in units of energy x flux rather than a ratio to a power law. The 
solid line is for the example shown in the top plot where r a — 60 
and the dotted line is for r = 10. When r = 60, the red wing of 
the iron line dominates during the fall meaning that the solid line 
in the bottom plot dips below zero for 5.4 > E > 6.4. During the 
fall, the blue wing dominates which gives rise to the hump in the 
6.4 > E > 7.4 region. Due to Doppler (and relativistic) boosting, 
the blue hump is larger than the red dip. When r — 10, the inner 
regions of the disc are moving much faster than the r = 60 case 
and therefore boosting is a much more significant effect. So much 
so, in fact, that the red wing of the iron line never dominates over 
the blue wing, even during the fall. The dotted line in the bottom 
plot therefore contains no red dip but only a blue hump. The peak 
of the blue hump is lower for r a = 10 than for r a — 60 but the 
area under the line is greater. This is because the iron line is more 
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heavily smeared in the r = 10 case, again due to faster orbital 
motion closer to the black hole. 

For both the LHS and the SIMS, the difference in iron line pro- 
file between the QPO rise and the QPO fall is significant, offering 
the possibility of direct observation for a range of spectral states. 
Note that this association of the rise with the bluest profile and the 
fall with the redest profile is robust as long as we are confident that 
the top (pole) of the flow is brighter than the sides. Because type-B 
QPOs provide a far cleaner signal than type-C QPOs, which are al- 
ways coincident with broad band variability, it will be easier to ob- 
serve this effect for a source in the SIMS. However, the QPO phase 
dependence of the iron line is particularly distinctive for the LHS 
model. An enhanced blue wing on the QPO rise, as predicted for 
the SIMS model, may feasibly be produced by some other process. 
A dominant red wing on the QPO fall and an enhanced blue wing 
on the rise, as predicted for the LHS model, can only realistically 
be produced by precession and a large truncation radius. Moreover, 
an observation showing that the difference spectrum changes be- 
tween states as we predict (i.e. the bottom plot of Figure [Tot would 
surely provide excellent evidence, not only of the precession model, 
but also that the truncation radius moves between the LHS and the 
SIMS. In the next section, we assess the likelihood of achieving 
such observational confirmation. 

5.2 Simulated observations 

We test the feasibility of observation directly by simulating phase 
resolved spectra using the FTOOL FAKEIT. This adds Poisson noise 
to a model before subtracting a representative background and de- 
convolving around a given response matrix. We simulate LHS spec- 
tra for 50 phase bins evenly spaced in precession phase angle, 7. We 
assume 100s exposure for each phase bin. This corresponds to 50 x 
100s = 5ks of good time. We sort the simulated data into four phase 
bins just as we did with the model. For the simulated data, there 
is just one QPO cycle with a long exposure but, for observational 
data there will be many short exposure QPO cycles to average over. 
As long as any fluctuations in the accretion geometry over this time 
are varying around an average value, the two processes should be 
equivalent to a good approximation. 

The top left plot in FigurefTTIshows the result of simulating the 
response of the Rossi x-ray timing explorer (RXTE) proportional 
counter array (PCA; top layer, detector 2). We unfold the spectrum 
around a flat power law and, as for the model, take the ratio to a 
power law with photon index T — 1.6. We use the same model as 
that shown in the top plot of Figure [Tol i.e. r = 60, <f>i = 90°, 
6i = 60°. Again, the green points are the trough, the blue points 
are the rise, the black points are the peak and the red points are 
the fall. Although a shift in line energy is visible between the rise 
and fall spectra, it is unlikely to be statistically significant due to a 
high noise level and low spectral resolution. The two observations 
of GRS 1915+105 studied by Miller & Homan (2005) were both 
observed with RXTE and, as such, the data were of a comparable 
quality to our simulation. They fit the QPO peak and trough spec- 
tra with a simple continuum model plus a Gaussian function for 
the iron line. When allowed to be free in the fits, the centroid en- 
ergy of the Gaussian was higher for the trough spectrum than for 
the peak spectrum in both observations. However, they were also 
able to achieve statistically acceptable results by fixing the centroid 
energy to the value measured for the total spectrum. Therefore, al- 
though there is some evidence that the line energy shifts, it is by 
no means statistically significant. It should be possible to achieve a 
slightly more significant result with RXTE data by comparing the 



rise and fall phases rather than the peak and trough, but this is al- 
ways marginal in practice due to the limited energy resolution of 
RXTE fast timing modes. 

The top right plot of Figurel 1 1 Ishows the same thing but for the 
XMM Newton European photon imaging pn camera (EPIC-pn). The 
Poisson noise level seems to be marginally worse compared with 
the simulated PCA data. Although the spectral resolution of the 
EPIC-pn is far better than that of the PCA, its effective area is less 
(~ 0.05m 2 compared with ~ 0.12m 2 ) meaning that we require a 
very heavy re-binning to get a reasonable signal to noise. Therefore, 
it may prove difficult to observe this effect using either RXTE or 
XMM Newton. However, the number of counts in the rise and fall 
phase bins could be maximised by halving the peak and trough 
phase bins and adding them to either the rise or the fall (i.e. the 
first half of the peak phase becomes part of the rise and the second 
half becomes part of the fall). 

A longer exposure is required to reduce the counting errors. 
In the bottom right hand panel of Figure [TT] we plot the result of 
assuming a lOOks exposure for the EPIC-pn. Encouragingly, we see 
that the dominant red wing in the falling phase is indeed resolved. 
However, over such a long exposure time, parameters such as r a 
may have systematically moved and so care must be taken to take 
this into consideration. 

The size of the effect is also dependent on our assumptions. A 
smaller flow scaleheight would increase the size of this effect be- 
cause the flux emitted from the poles of the flow would be an even 
greater fraction of the flux emitted from the entire flow. Frame drag- 
ging could therefore have a larger effect on the iron line than we 
predict here making it easier to observe with current instruments 
than our simulations imply. However, it also must be noted that 
the continuum will be more complicated than we assume here with 
some QPO phase dependent spectral pivoting resulting from a vari- 
ation in the flux of disc photons incident on the flow. This will make 
observation harder. 

The bottom left hand plot of FigurefTTIshows the potential im- 
pact of the proposed mission LOFT (the large observatory for x-ray 
timing). We use the 'required' response of the large area detector 
(LAD), which is the principle instrument of the mission. Because 
the LAD has an exceptionally large effective area (10- 12m 2 ), the 
results are far clearer than those provided by current missions. In 
fact, the noise level is so low with LOFT, it would be possible to 
constrain spectra for far more than four phase bins. We could also 
constrain these spectra for less than 5ks good time, meaning that 
we could conduct detailed studies of the evolution of the phase re- 
solved spectra. 

5.3 RMS spectum 

Since we calculate 50 spectra for both the LHS and SIMS models, 
it is simple to calculate the rms spectrum of the QPO. This is sim- 
ply the standard deviation of each energy channel in absolute units 
(i.e. not divided through by the average). Figure [T2] shows this for 
the LHS model (top) and the SIMS model (bottom) with the mean 
spectrum plotted in black and the QPO spectrum plotted in red. 
Since the QPO spectrum is fairly sensitive to model assumptions, 
it provides a good way to constrain model parameters against ob- 
servation. For the models we use here, the misalignment angle f3 is 
large and thus we see reflection features in the LHS QPO spectrum 
as the amount of reflection changes with QPO phase. 

By contrast, in the SIMS, the extent of the flow is so small 
(ri — 7 and r = 10) that even this large misalignment angle does 
not give rise to significant variability in the total reflection fraction. 
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Figure 11. Simulated observations of the phase binned spectra shown in Figure [TOl with r = 60, (/>{ = 90° and 9i = 60°. These spectra are unfolded around 
a flat power law and plotted as the ratio to a power law with T = 1.6 and unity normalisation. Again the four phase bins are for the QPO minimum (green), 
rise (blue), maximum (black) and fall (red). Observed with the RXTE PCA or the XMM Newton EPIC-pn, it is difficult to see by eye the difference in iron 
line peak energy between different phase bins. In contrast, the Athena WFI recovers the model well and the LOFT LAD does so with an exceptionally high 
precision. 



Previous rms spectral analysis of the QPO have not looked at this 
in detail (e.g. Sobolewska & Zycki 2006). We plan to address this 
issue in a future work (Axelsson et al in preparation). 



6 CONCLUSIONS 

The truncated disc / precessing inner flow model for the spectral 
timing properties of XRBs predicts a QPO phase dependence of 
the iron line profile. This results from the inner flow preferentially 
illuminating different regions of the disc as it precesses. When the 
brightest region of the disc is moving towards us, the iron line will 
be blue shifted and boosted. When the brightest region is receding, 
the iron line will be red shifted. As the illumination pattern rotates 
around the disc, the iron line rocks between blue and red shift. This 
process always happens in a particular order with the most heavily 
blue shifted iron line profile following the QPO trough and the most 
heavily red shifted iron line profile following the QPO peak. It is 



possible to isolate the peaks and troughs in a light curve using a 
simple flux selection. The rising phase, which follows the trough, 
is predicted to have the bluest iron line and the falling phase, which 
follows the peak, is predicted to have the reddest iron line. 

We predict this QPO phase dependence of the iron line pro- 
file to be present for a large range of spectral states (and therefore 
truncation radii). This means that it may be best to search for the ef- 
fect in spectra containing type-B QPOs which have very little broad 
band variability associated with them and therefore provide a much 
cleaner signal than type-C QPOs. However, the nature of the iron 
line phase dependence changes with truncation radius. When it is 
large, the red wing can dominate over the blue wing during the fall 
from QPO peak to trough. When it is small, Doppler and relativistic 
boosting from the rapidly moving inner regions of the disc means 
that the red wing can never dominate over the blue wing. The char- 
acteristic shape of the difference spectrum between rise and fall 
should therefore change as the spectrum evolves from the LHS to 
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Figure 12. Mean and QPO spectra for the LHS (top) and SIMS (bottom) 
models. The QPO spectum is calculated by measuring the standard devi- 
ation of each energy channel around the mean value across 50 values of 
precession angle. 



the SIMS. The dominant red wing of the QPO fall spectrum in the 
LHS (the 'red dip' in the difference spectrum) is the most unique 
model prediction but if we wish to observe this, we must disen- 
tangle the underlying QPO signal from the broad band noise. This 
will be the subject of a future paper. An observation of the effect 
in both states, along with confirmation that the difference spectrum 
changes with state, would constitute excellent evidence, not only 
of the precession model, but also that the truncation radius moves 
between the LHS and the SIMS. 

Quasi-periodic shifting of the iron line peak energy is a unique 
prediction of the Lense-Thirring precession model for the low fre- 
quency QPO in XRBs. We have shown that it may be possible to 
observe such an effect with current missions, but that LOFT will be 
able to measure this with precision, enabling us to place accurate 
constraints on the accretion geometry. 
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APPENDIX A: GEOMETRY 

In order to perform our calculations, we must define some vectors 
using the coordinate system outlined in Figure[T] We represent the 
x, y and z axes with the standard i, j and k unit basis vectors. It 
then follows from FigureQjthat 



J 



BS 



— sin P j + cos P k 
cos (3 j + sin /3 k. 



(Al) 



The three vectors i, e and J_ BS therefore form a right handed Carte- 
sian coordinate system: the disc basis vectors. We can define a vec- 
tor, Td r d , which points from the origin (the black hole) to any point 
on the disc where 



: cos (j>di J r sin <p& e. 



(A2) 



Note, because the disc is razor thin, there is no j_ BS component (i.e. 
J_ BS .£ d — 0) and <j>d is simply the angle between r d and the x-axis. 
We also define a vector pointing from the origin to the observer 
using the disc basis vectors 



S_ = sin 0i cos 4>i i + sin Qi sin 4>i e + cos Qi j_ BS . 



(A3) 



In order to describe points on the surface of the flow, we must 
define flow basis vectors. The 'z-axis' of this right handed coor- 
dinate system is J_f low which precesses around k as illustrated in 
Figure [T] The other two basis vectors, x^ and $ , must therefore 
also precess with the flow. We use 



Because our flow is elliptical with semi-minor axis in the 
j-fiow direction and semi-major axis in the a = cos</>/ x_j + 
sin cf>f y direction, the distance from the origin to any point on 
the surface is 



r h 



^(h sine?/) 2 + (r D cost?/) 2 



(A6) 



Because r/ is uniquely determined by Of, we can define dr = 
\r(6 f )-r(e f +d0f)\. 

We need to be able to write down the unit vector normal to the 
flow surface. We can do this using a few identities. Imagine a trian- 
gle drawn between the two focuses of the ellipse, Fi and F%, and 
any point on the circumference of the ellipse, P. We know that the 
distance from the origin to either focus is / = \/r 2 — h 2 and also 
that the three sides of the triangle add up to 2r + 2f. We can define 
the angle between the line from P to F\ (P F± ) and the line from 
P to F2 (P F2) as ip. We know that the surface area unit vector, A, 
goes directly between these two lines such that the angle between 
—A and each line is i[i/2. We can say that A points from some point 
x a to the point on the flow surface, P, in such a way that this con- 
dition is satisfied. Say that d is the distance from P to F2 and Q is 
the angle between the lines F2 Fi and F2 P. We can use the cosine 
rule a few times to show that d - 



f 2 + r'j — 2/r/ sin 8f and 
cos Q, = (f 2 — r 2 + r d) /(/ d). It is then possible to show that 



A = 



T 'f Lf — X oO L 



^jx 2 + r 2 — 2x rf sin Of 



where 

cosi/> = 
and 



2r 2 + d 2 - 2r d - 2f 2 
d(2r a - d) 

dsin(V>/2) 
~ sin(7r-V>/2-n)' 



(A7) 



(A8) 



(A9) 



We will also need to define a vector which points from a given 
point on the flow to a given point on this disc. This can be written 

as 

C£ = -rftf + rdtt- ( A1 °) 

From this, it is simple to show that the distance between the two 
points is 



f2 2 . 2 

C = Tf + r d - r f r d r f .r d r d . 



(All) 



All of these vectors will become very useful for the following sec- 
tions. 



J 



flow 



= cos 71 + sin 7 j 

= — cos P sin 7 i + cos P cos 7 j + sin P k 
= sin P sin 71- sin P cos 7 j + cos P k, 



(A4) 



such that Xf — i when 7 = but, as the precession angle unwinds, 
the axes move. We can then specify a point in the flow with the 
vector rf r_f where 

f_f — sinf?/ cos 4>f x_j + sin 6*/ sin0/ y + cos 8f J_f low - (A5) 

Here, Of is the angle between and J_f low and cj>f is the angle 
between r,(0f — n/2) andi. 



APPENDIX B: DISC IRRADIATION CALCULATIONS 

So, we need to calculate what luminosity a disc element with sur- 
face area dAd = rddcfradrd will intercept from a flow surface ele- 
ment emitting a luminosity dL over a semi-sphere (because it only 
emits away from the rest of the flow). We can then integrate over 
all flow elements to work out the total flow luminosity that the disc 
element intercepts. For the disc patch to see anything at all from a 
given flow element, it must pass two tests. First, does it lie in the 
unit semi-sphere of the flow element; i.e. is A.C, > 0. Also, because 
we are viewing the top of the system (9i ^ 90°), we only see lu- 
minosity which has reflected off the top of the disc. Therefore, we 
only count luminosity incident on the top of the disc in our integral. 
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This means we require C-<Lbs < 0. If one of these conditions isn't 
met, the luminosity intercepted by the disc element is dL r — 0. If 
both are, we have 



dL r = - — — 1 



y )dA c 



2ttC 2 



dh. 



(Bl) 



We see that, the amount of luminosity intercepted depends on the 
projected area of the disc patch as seen by the flow element. If 
the patch is face-on as seen by the flow, C-J-bs ~ 1 an d me P ro " 
jected area is dAd- This area reduces as the patch turns away from 
the emitting flow element. The total luminosity incident on a disc 
patch is calculated by adding up the contribution from every flow 
element. 



APPENDIX C: FLOW MODULATION CALCULATIONS 

We now need to calculate how much luminosity a telescope with 
effective area A e f / will intercept from a given flow element in or- 
der to again integrate over the whole flow. For the telescope to see 
any luminosity at all, two tests must again be passed. First of all, 
the viewer must be in the unit semi-sphere of the flow element. This 
means we require A.S_ > 0. We also wont see anything if the emis- 
sion is blocked by the disc. We know the emission definitely wont 
be blocked by the disc if the flow element is above the disc; i.e. 
Lf-J-BS -* 0- Even if the element is below the disc plane, we still 
might be able to see through the hole in the centre of the disc. So, 
imagine a point on the flow which is below the disc plane, emitting 
along the vector S_. At some point it will intercept the disc plane. 
The distance between the flow element and the point where the vec- 
tor crosses the disc plane is f . This point will be a distance rd from 
the origin. We can write 



(,§.= -r f r f +r d r d . 

Dotting both sides with J_ BS and rearranging gives 



c = 



r f Lf-J-BS 



cos 0i 
We then know that 



2 

rd 



C+rj+ 2Cr f S.r f 



(CI) 



(C2) 



(C3) 



If r d < Tg, we still see the flow element through the hole in the 
disc. If not, it is hidden by the disc. 

So, if the unit-sphere and disc obstruction tests are not passed, 
the luminosity intercepted by the telescope is dL t s = 0. Other- 
wise, this is 



dL bs — 



i eff 



dL 



(C4) 



2ttL> 2 ' 

where D is the distance to the source. Note, because the telescope 
is so far away and is pointed straight at the black hole, we can say 
that the projected area of the telescope as seen by any flow element 
is A e ff. We then just set A e ff/(2nD 2 ) = 1, because it only tells 
us about normalisation, and sum up the contribution from each flow 
element. 



points of the disc will therefore travel towards the observer at a ve- 
locity of ±Vk sin Si. This means that a photon emitted with energy 
Eem with be red shifted by 



E em /E obs = (l-3/r d )- 1/2 



Ml + tan 2 Co)-2] 1 /2 



,(D1) 



where 



+ cos 2 d> sin 2 t 



-1/2 



tan£ D = cos cj> sin#i(l — cos 2 4>sin 2 9i) 



-1/2 



(D2) 



(Fabian et al 1989; 2000). 

For a given precession angle, 7, the flow luminosity incident 
on a disc patch described by rd and <j>d is L r (rd, <j>d)- If this lumi- 
nosity were all emitted at energy E em , the observer would see a 
luminosity, all at E b s , of 



dL obs w L r (r d ,cj> d )(E oba /E e m) cosdi. 



(D3) 



Here, the approximations come from assuming light to travel in a 
straight line. Throughout this paper, we ignore gravitational light 
bending thus taking these to be good approximations. This should 
be apropriate since the inner radius of the flow is assumed to be 
Ti = 7 throughout and light bending effects outside of this radius 
will be minimal. The total observed luminosity as a function of 
energy is calculated by summing the contribution from each disc 
patch. As the flow precesses and the function L r (rd, 4>d) evolves, 
the observed iron line profile will change. 



APPENDIX D: IRON LINE PROFILE CALCULATIONS 

A disc element at r d is rotating with Keplerian velocity Vk. An 
observer at 8i, 4>i then sees the disc patch travelling towards them 
at a velocity of v = Vk sin cj> sin 9i where <f> = <f>i — <f>d- The tangent 
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